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A hybrid simulation code is developed that is suited for fast one-dimensional simulations of shower profiles, 
including fluctuations. It combines the Monte Carlo simulation of high energy interactions with a fast numerical 
solution of cascade equations for the resulting distributions of secondary particles. Results obtained with this 
new code, called conex, are presented and compared to corsika predictions. 



1. INTRODUCTION 

Ultra-high energy cosmic rays (UHECR) are 
investigated by measuring the characteristics of 
the secondary particle cascades that they induce 
in the Earth's atmosphere. The information ob- 
tained from these extensive air showers (EAS) is 
used to infer the properties of the primary parti- 
cle, relying on a proper theoretical description of 
the cascade processes. 

The most natural way to get detailed informa- 
tion on the atmospheric particle cascading seems 
to be a direct Monte Carlo (MC) simulation of 
the EAS development, like it is done, for exam- 
ple, in the CORSIKA program pp. Nevertheless, 
for primary particles of very high energies, this 
is not a viable option because of unreasonably 
large calculation time required. The situation can 
be improved by applying some weighted sampling 
algorithms, like the so-called "thinning" method 
[21 . Although this approach allows the reduction 
of EAS calculation times to practically affordable 
values, it has limitations. The summation of par- 
ticle contributions with very large weights creates 



*speaker, e-mail: Tanguy.Pierog@ik.fzk.de 
tNow at ° 



significant artificial fluctuations for EAS charac- 
teristics of interest (3| . 

A possible alternative procedure is to describe 
EAS development numerically, based on the so- 
lution of the corresponding cascade equations 
|4I5I6) . Combining this with an explicit MC sim- 
ulation of the most energetic part of an EAS al- 
lows one to obtain accurate results both for aver- 
age EAS characteristics and for their fluctuations 
Uj. In this article we report on the development 
of an extensive air shower calculation program 
of such a type: CONEX. The MC treatment of 
above-threshold particle cascading is treated in 
the standard way and does not differ significantly 
from e.g. CORSIKA procedure ^J. On the other 
hand, the numerical description of lower energy 
sub-cascades is based on the solution of hadronic 
cascade equations, using an updated algorithm of 
Ref. [S] and a newly developed procedure for solv- 
ing electro-magnetic (e/m) cascade equations. 

The outline of the paper is as follows. Sec- 
tion 12 describes the calculation scheme and its 
basic procedures. In Section |31 we present exam- 
ples of EAS characteristics obtained with CONEX 
and investigate the accuracy of the predictions 
comparing the hybrid approach with pure MC or 
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numerical procedures, and with the results of the 
CORSIKA code. Finally, Section^contains a sum- 
mary of our work and discusses both potential 
applications of the program and the prospects for 
its further development. 

2. SCHEME 

A hybrid air-shower calculation scheme consists 
of two main stages: an explicit MC simulation of 
particle cascade at energies above some chosen 
threshold Ethi (typically a factor of 100 smaller 
than the energy of the primary particle) and a so- 
lution of nuclear-electro-magnetic cascade equa- 
tions for sub-cascades of smaller energies. Both 
MC and numerical parts model the same physics, 
but the numerical solution to the cascade equa- 
tions is calculated only along the direction of the 
shower axis. 

In the hadronic cascade one follows the prop- 
agation, interactions and decays (if relevant) of 
(anti-)nucleons, charged pions, charged and neu- 
tral kaons. All other types of hadrons produced 
in the interactions and decays are assumed to 
decay at the place. Particle interactions in the 
MC part are treated within a chosen high energy 
hadronic interaction model (optionally nexus |S] 
or QGSJET "S") and decays are simulated using the 
corresponding routines of the nexus model. The 
same models are applied to pre-calculate average 
secondary particle spectra for later use by the nu- 
merical scheme. Optionally, below some energy 
100 GeV one employs the GHEISHA model 
|10l . Ionization energy loss of charged hadrons is 
approximated by the Bethe-Bloch equation. 

The MC treatment of the e/m cascade is re- 
alized by means of the EGS4 code |T1|, sup- 
plemented by an account for the Landau- 
Pomeranchuk-Migdal (LPM) effect for very en- 
ergetic electrons (positrons) and photons. The 
numerical part is based on the same interaction 
processes as the MC one, using Bethe-Heitler 
cross sections for bremsstrahlung and pair pro- 
duction with corrections at low energy according 
to Storm and Israel, the Klein-Nishina formula 
for the Compton process, accounting for M0ller 
and Bhabha processes as well as for positron- 
electron annihilation (see, for example |12lll| ). 



Both LPM-effect and photo-effect are neglected in 
the numerical part as it is not supposed to be used 
in the energy ranges where the corresponding pro- 
cesses give an essential contribution. Ionization 
losses of electrons and positrons are described by 
the Bethe-Bloch formula with corrections due to 
the density effect. 

In general, an individual shower is simulated 
as follows. One starts with the primary particle 
of given energy, direction, at a given initial posi- 
tion in the atmosphere. These define the shower 
axis along which the shower will be calculated 
using slant depth. For a primary hadron initi- 
ating a shower one simulates the hadronic cas- 
cade explicitly until all produced hadrons fall be- 
low the energy threshold E'thr- The characteris- 
tics of all sub-threshold hadrons and e/m par- 
ticles (type, energy, and slant depth position) 
are written to corresponding stacks to form the 
"source terms" for the numerical solution of the 
cascade equations. The above-threshold e/m par- 
ticles are transferred to the EGS4 program, where 
the simulation of the e/m particle cascade is per- 
formed in a similar way, with all sub-threshold 
e/m particles being added to the e/m source term. 
As the next step, the hadronic cascade at ener- 
gies below Ethr is described numerically based on 
the solution of corresponding cascade equations 
and the initial condition specified by the source 
term. The results are discretized energy spectra 
of hadrons of different types at various depth po- 
sitions. All sub-threshold e/m particles produced 
at this stage are added to the e/m source term. 
Finally, sub-threshold e/m cascades are described 
numerically based on the solution of correspond- 
ing e/m cascade equations with the initial con- 
ditions set by the corresponding source term. In 
case of the primary particle being a photon or an 
electron the simulation process consists just from 
e/m MC cascading for above-threshold particles 
and numerical treatment of the process for sub- 
threshold particles. No feed-back from the e/m 
to the hadronic part is considered. 

3. RESULTS 

Currently CONEX can be used for the simula- 
tion of showers initiated by any primary particle 
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Figure 1. Average longitudinal profiles of charged 
particles of energies above 1 MeV for proton- 
initiated vertical showers of energy Eq = 10^^ eV 
as calculated both in the hybrid scheme (dashed 
line) and using pure MC (points) or numerical 
(CE) approaches (full line), compared to COR- 
SIKA simulation (stars). 
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Figure 2. Particle energy spectra of electrons and 
positrons for the depths X — 700 and 1000 g/cm^ 
for lO^'* eV proton-initiated vertical showers of 
energy as calculated both in the hybrid scheme 
(dashed line) and using pure MC (points) or nu- 
merical (CE) approaches (full line), compared to 
CORSIKA simulation (stars). 



and energy accepted by the high energy inter- 
action model implemented (7, proton, iron and 
more up to 10^^ eV). The basic output consists 
of the longitudinal profile along the shower axis 
for all produced particles (e/m and hadronic) ex- 
cept for the muons (under development). Ad- 
ditionally, Gaisser-Hillas fit parameters, longitu- 
dinal energy deposit profiles and energy spectra 
for 3 different slant depth levels can be provided. 
The geometry of the model allows the simulation 
of showers at any angle up to 90° assuming US 
standard atmosphere, but can easily be adapted 
to any atmosphere and upward going showers. On 
the other hand, since numerical calculations are 
only done along the shower axis, no lateral distri- 
bution is available. 

In Fig. n the average longitudinal profile of 
charged particles for 10^^ eV proton-initiated ver- 
tical showers as calculated both in the hybrid 
scheme (i?thr = 10'' GcV) and using pure MC 
{Ethr = 1 GeV) or numerical (£:thr = 10^ GeV) 
approaches with GHEISHA and QGSJET is com- 
pared to CORSIKA simulations. All methods are 
fully compatible with CORSIKA results. Fig. El 
shows in more detail how well the different meth- 
ods and the CORSIKA results agree in the de- 
scription of the energy spectra of e/m particles 



around the maximum {X = 700 g/cm^) and at 
the ground level {X = 1000 g/cm^) for 10^* eV 
proton-initiated vertical showers. 

The most important feature for the hybrid 
scheme is to reproduce the event- by- event fluc- 
tuations of the shower development as a function 
of the primary mass. In Fig. 13 fluctuations of the 
shower maximum depth Xmax around the mean 
shower maximum depth (ATmax} for a primary 
energy of lO^^eV are shown. Proton and iron- 
initiated showers were simulated with CONEX and 
are compared with CORSIKA results. A perfect 
agreement is observed for both primary particles. 

In Fig. 01 the average depth of maximum, 
(ATmax) , of proton and iron-induced air showers is 
shown as function of the energy for various mod- 
els and data. The the predictions of the CONEX 
and CORSIKA simulation codes agree within the 
statistical uncertainties. 

4. SUMMARY 

CONEX is a new, reliable tool for fast one- 
dimensional air shower simulations. Realistic 
shower profiles, including fluctuations that lead, 
for example, to showers with two maxima, can be 
generated for particles of primary mass, energy 
and incident angle as relevant to UHECR studies. 
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Figure 3. Fluctuations of the shower maxi- 
mum depth Xniax around the mean shower max- 
imum depth (Xniax) for 10^^ eV proton and iron- 
initiated showers simulated with CONEX (lines) 
and compared with CORSIKA results (points). 



Figure 4. (^max) of proton and iron-induced 
showers for different models together with data. 
CORSIKA shower simulations are based on 500 
showers for each energy point. 



With a calculation time of about 2 minutes per 
shower (depending weakly on the primary energy 
and mass), CONEX provides necessary informa- 
tion for applications such as fluorescence detector 
simulations or theoretical studies of shower pro- 
files and Xinax distributions. 

Work is in progress to include muon calcula- 
tions and new high and low-energy hadronic in- 
teraction models in CONEX. Furthermore it is 
planned to link CONEX to CORSIKA to provide full 
3D information and low energy particle tracking 
for surface detector analysis or shower radio emis- 
sion generation. 
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